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ABSTRACT 

We present the first interferometric polarization map of the W3(0H) massive 
star-forming region observed with the Submillimeter Array (SMA) at 878 /im 
with an angular resolution of l'/5 (about 3 x 10 3 AU). Polarization is detected 
in the W3(H20) hot core, an extended emission structure in the north-west of 
W3(H 2 0), and part of the W3(0H) ultracompact Hn region. The W3(H 2 0) 
hot core is known to be associated with a synchrotron jet along the east-west 
direction. In this core, the inferred magnetic field orientation is well aligned with 
the synchrotron jet and close to the plane of sky. Using the Chandrasekhar- Fermi 
method with the observed dispersion in polarization angle, we estimate a plane- 
of-sky magnetic field strength of 17.0 mG. Combined with water maser Zeeman 
measurements, the total magnetic field strength is estimated to be 17.1 mG, 
comparable to the field strength estimated from the synchrotron model. The 
magnetic field energy dominates over turbulence in this core. In addition, the 
depolarization effect is discerned in both SMA and JCMT measurements. Despite 
the great difference in angular resolutions and map extents, the polarization 
percentage shows a similar power-law dependence with the beam averaged column 
density. We suggest that the column density may be an important factor to 
consider when interpreting the depolarization effect. 

Subject headings: ISM: individual (W3(H 2 0)) — ISM: magnetic fields — Tech- 
niques: polarimetric — stars: formation 



1. Introduction 

One of the outstanding problems in star formation concerns the influence of magnetic 
fields to the dynamics and timescales (Shu et al. 1999; McKee & Ostriker 2007). The rela- 
tive importance between magnetic fields and turbulence affects cloud evolution and cluster 
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formation (McKee & Ostriker 2007). Ambipolar diffusion has been proposed to regulate 
core collapse (Shu, Adams, & Lizano 1987). In accretion flows, magnetic effects such as 
magnetic tension, braking, and reconnection, or non-ideal magnetohydro dynamic (MHD) 
effects, can drive the evolution of the formation process (Shu, Li, & Allen 2004; Li, McKee, 
& Klein 2011). The residual magnetic fields in the circumstellar environs further assist to 
collimate supersonic outflows (Shu et al. 2000; Shang, Li, & Hirano 2007). Intringingly, 
radio emission with negative spectral index has been reported in a few protostellar jets, e.g. 
HH80-81 (Marti, Rodriguez, & Reipurth 1993; Carrasco-Gonzalez et al. 2010), Cepheus A 
(Garay et al. 1996), W3(H 2 0) (Reid et al. 1995; Wilner, Reid, & Menten 1999). Such non- 
thermal emission resembles jets driven by active galactic nuclei and suggests a connection 
with synchrotron jets. 

It has become increasingly clear that dust grains are aligned by magnetic fields not 
only in the diffuse interstellar medium but also in dense cloud cores or circumstellar disks 
(Myers & Goodman 1991; Ward-Thompson et al. 2000; Lazarian, Goodman, & Myers 1997; 
Matthews et al. 2009; Dotson et al. 2000, 2010). Although grain alignment mechanisms 
may be less assured, polarization remains one of the most informative techniques to trace 
magnetic fields (Lazarian 2007). Complimentary to polarization from selective extinction of 
background starlight, polarized emission arising from magnetically aligned dust grains in the 
millimeter or submillimeter wavebands is often used to map the two-dimensional morphology 
of the three-dimensional (3D) magnetic fields projected on the plane of sky, B ± . Depending 
on the type of polarization, the observed orientation is supposed to be either parallel to B± if 
observed in dust absorption against background sources or perpendicular to Bj_ if observed 
in dust thermal emission (see the review by Hildebrand 1988). 

The bright ultracompact (UC) Hll region, W3(OH), marks a nearby (d = 2.04 kpc; 
Hachisuka et al. 2006), well-studied massive star- forming region in the course of developing 
an OB stellar group. At least seven H II regions at different evolutionary stages were found 
within 30" of W3(OH) (Harten 1976). The most recent star formation takes place in a lumi- 
nous (L ~ 1O 4 L ) hot core, W3(H20) (a.k.a the Turner- Welch object), which is associated 
with water masers and appears as a warm, dense molecular clump about 6" east to W3(OH) 
(Turner & Welch 1984). Millimeter and submillimeter interferometric studies resolved the 
dust emission into two components, A and C, and suggest an accreting proto-binary system 
(Wyrowski et al. 1999; Chen et al. 2006; Zapata et al. 2011). Proper motions of water 
masers trace a bipolar outflow centered at W3(H 2 0)A (Alcolea et al. 1993; Hachisuka et al. 
2006), which also drives a synchrotron jet along the east- west direction (Reid et al. 1995; 
Wilner et al. 1999). Here we perform interferometric polarization observations toward the 
W3(OH) massive star-forming region to map the morphology of the projected magnetic fields 
in this region. Our study is enabled by the Submillimeter Array 4 (SMA) with its polarimeter 
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system (Marrone et al. 2006; Marrone & Rao 2008). 

2. Observations and Data Reduction 

Observations using quarter-wave plates installed on the SMA were carried out with 
seven antennas in the extended configuration on 2007 August 10 and with six antennas in the 
compact configuration on 2008 September 20. The phase tracking center is at (a, <5)(J2000) = 
(2 h 27 m 3.870 s , +61°52'24.60"). The system temperature varied from 150 to 580 K during the 
extended array observation and from 110 to 400 K during the compact array observation. The 
correlator was set to have a spectral resolution of 0.8125 MHz, corresponding to 0.7 kms -1 , 
and 2 GHz bandwidth in each sideband. The projected baselines have a range of 8 — 207 kA, 
which is insensitive to structures larger than 11" (Wilner & Welch 1994). The full- width at 
half power (FWHP) of the primary beam is roughly 37". 

Data inspection, calibration, and imaging were performed with the MIRIAD package. 
The instrumental polarization response, i.e. "leakage," for each antenna was calibrated 
by observing strong quasars, 3clll and 3c454.3, over a wide hour angle range for a good 
parallactic angle coverage. The bandpass calibration was performed with 3c454.3. Ganymede 
and Callisto were used to set the flux density scale, which is estimated to be accurate 
within 15%. Regularly interleaved observations of the nearby quasar 3c84 were conducted 
to monitor the atmospheric variations. Atmospheric complex gains were first solved with 
respect to 3c84 to make initial Stokes / images for individual tracks. The gain solutions 
were further improved by one iteration of phase-only self-calibration using the initial Stokes 
I images as models. To image continuum emission, care has been taken to remove channels 
with line emission. Both sidebands were used to generate line-free continuum maps in a 
multi- frequency synthesis with an effective central frequency of 341.5 GHz. 

To obtain the best signal-to- noise ratio, maps of Stokes / (total intensity), and Stokes 
Q and U (linear polarization) were made with visibilities weighted by the associated sys- 
tem temperature using natural weighting. This results in a synthesized beam of I'M x 1"5 
with position angle (P.A.) of —14°. For easy comprehension and convenience of subsequent 
Nyquist sampling, we smooth all the maps with a circular Gaussian beam of 1"5. The re- 
sultant maps have an rms noise level of 07 = 27 mJybeam -1 for the Stokes / map and 
&qu — 3 mJybeaur 1 for maps of Stokes Q and U. The dynamic range of the Stokes / map 
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is limited at a maximum intensity value of 120 07 due to the u-v coverage and calibration 
accuracy. Once Stokes /, Q, and U are known, the polarization percentage, p, and the po- 
larization position angle, 9, can be deduced (Matthews & Wilson 2000). The raw polarized 
intensity, I P)Tam = (Q 2 + f/ 2 ) 1 / 2 , is necessarily positive and biased by the rms noise of Stokes 
Q and U, uqu- Hence, a debiased polarized intensity is used, I p = (Ip >raw — Vqu) 1 ^ 2 (Leahy 
1989; Vaillancourt 2006) when computing the polarization percentage, p = I p /I. The po- 
larization position angle is determined by 6* b s = (1/2) tan _1 (?7/Q), with an uncertainty of 
cfyobs = (1 / 2) (o-Qu/Ip) (180° /it). Note that in linear polarization measurements, values of 
#obs different by 180° are identical. After generating maps of Stokes I, Q, and U, we image 
the distribution of polarization, including I p , p, and 9 ^ s , for regions where l v > 3ctqu and 
I > 3ctj. Tabled] lists the Nyquist sampling of the polarization map with a grid step of ft' 75. 



3. Results and Discussion 

The 878 /xm polarization map is shown in Figure [TJ Line segments indicate the orienta- 
tion of polarization angle, # b s , with their lengths proportional to the polarization percentage, 
p. Both the W3(H20) hot core and the W3(OH) UC Hn region are clearly detected in the 
Stokes I map. Polarization is detected in W3(H 2 0), W3(OH), and their vicinity, particu- 
larly in a faint, extended component in the north-west of W3(H 2 0). In the W3(OH) UC Hn 
region, the total intensity is dominated by the free-free continuum, which is supposed to be 
unpolarized. Additional care is needed to estimate the relative contributions from dust and 
free-free emission to the total intensity in W3(OH), and this will be discussed elsewhere. 



3.1. Morphology and Properties of Detected Polarization 

In Tabled], we group the detected polarization values based on their apparent association 
with distinct polarization subregions: W3(H 2 0), W3(H 2 0)-NW referring to the extended 
structure in the north-west of W3(H 2 0), W3(H 2 0)-E and W3(H 2 0)-W referring to the two 
small regions in the east and west of W3(H 2 0), respectively. With the exception of W3(H 2 0)- 
NW, there is a general east-west orientation of B± across the whole region (Figure [2]). 

In W3(H 2 0), the polarization orientation is remarkably uniform, 9 D b s = —12° — 20°. 
The mean polarization angle weighted by observational uncertainty is (0 o bs) = 0° ± 2°. We 
compute the standard deviation, 59 \, s , with respect to (0 o bs) and obtain a value of 11.4° 
(Table [2]). This observed dispersion, 89 Q b s , includes the measurement uncertainty, which is 
the root mean square (rms), Scrg^, of the observed uncertainties, ctyobs- Hence, the intrinsic 
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dispersion of polarization angle is given by 59 = (59 2 ohs - <Wobs) 1/2 = 9 - 4 ° and wil1 be used 
to estimate the magnetic field strength of B± (Sect. 13. 2j) . In addition, the distribution of 
polarization angle in W3(H 2 0)-NW (Figured]) also appears fairly uniform. We compute for 
W3(H 2 0)-NW a mean polarization angle of (0 o bs) = 54° ± 2° and an intrinsic polarization 
angle dispersion of 58 = 10.2°. 



3.2. The Magnetized Environment of the Synchrotron Jet in W3(H 2 0) 

The synchrotron jet exhibits a double-sided morphology along the east-west direction 
centered at the millimeter continuum peak W3(H 2 0)A in the VLA 3.6 cm map (Figure [2J 
Wilner et al. 1999). Our uncertainty weighted mean polarization angle of (^obs) = 0° ± 
2° implies an east-west orientation of B±, in good agreement with the direction of the 
synchrotron jet. In early VLA studies of the synchrotron jet (Reid et al. 1995; Wilner 
et al. 1999), an inhomogeneous model of a biconical geometry with opening angle of 0.1 rad 
successfully reproduced the observed properties. The model gave a total magnetic field 
strength of B(r) = 10 mG (r/0'.'2)~ a8 , where r is the angular distance from the center. 

The field strength of the projected magnetic field component, B±, can be estimated if 
the dispersion of the polarization angle, 59, the mass density, p, and the velocity dispersion 
along line of sight, 5v\ os , are measured. Considering a uniform magnetic field, Chandrasekhar 
& Fermi (1953; hereafter CF) first proposed a "polarization-dispersion" method to estimate 
B± based on transverse motions of Alfven waves, which are observable as dispersion in 
polarization angles, 58, from the mean orientation of the projected fields. If the velocity 
dispersion is isotropic and the transverse velocity dispersion is equal to 5vi os , one can obtain 
the plane-of-sky magnetic field strength 

= M^(±)(JL)»(-^)»(«!!- )( S 1Y\ (1) 

V0.5yVl.36/ Vl0 7 cm- 3 / V lkms /V 10 / 

where Q = 1 corresponds to the original CF formula, and \x is the mean molecular weight. In 
a later study, Ostriker, Stone, & Gammie (2001) performed 3D numerical MHD simulations 
with conditions in turbulent clouds and suggested a factor of Q ~ 0.5 applicable at varying 
observer orientations as long as 56 < 25°. Chen et al. (2006) have studied the physical 
conditions in W3(H 2 0) by analyzing the 1.4 mm and 2.8 mm continuum and the CH 3 CN 
lines. Since our polarization observations do not resolve individual protostellar envelopes, 
we simply take their mean value of wh 2 = 1.5 x 10 7 cm -3 and 5v\ os = 1.9 kms -1 based 
on the CH 3 CN spectra. Using Equation (pQ), we estimate the magnetic field strength to be 
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B± = 17.0 mG. Since 5# bs can be lowered by the finite angular resolution and averaging 
along line of sight, this estimate of B± should be considered as an upper limit. 

When Zeeman measurements are available to provide the line-of-sight magnetic field 
strength, B\\, it is possible to evaluate the total field strength with B = (B\ + B 2 ) 1 ^ 2 . This 
offers an independent method to estimate the total field strength in comparison with the 
synchrotron model. The Zeeman measurements in 22 GHz H2O masers in W3(H20) yield 
a postshock field strength of -B||, p0 st = +42.1 mG (Sarma et al. 2002). Depending on the 
degree of shock compression, the field could be amplified by a factor of about 20 (Sarma et al. 
2002) and renders an estimate for the preshock field strength of B\\ = +2.1 mG. Hence, an 
estimate of the total field strength is about B ~ 17.1 mG, with orientation close to the plane 
of sky. Though slightly higher, this total field strength is comparable to the value given by 
the synchrotron model (Reid et al. 1995; Wilner et al. 1999). 

Following Crutcher (1999), we estimate the observed-to-critical mass-to-magnetic flux 
ratio, (M/$ B ) obs/crit = (iV H2 /10 24 cm.- 2 )/(B/10 mG) = 0.54, for W3(H 2 0) with Nn 2 = 
9.2 x 10 23 cm -2 (Chen et al. 2006). This implies the core is magnetically subcritical and 
can be supported by the static magnetic field against the gravitational collapse. However, 
ongoing collapse in this region has been suggested by infall asymmetry in spectral profiles 
(Wu & Evans 2003; Wu et al. 2010). Since interferometric observations tend to under- 
estimate column density and B± is an upper limit, our estimate of (M/$s) obs / crit should 
be treated as a lower limit. Using the total field strength, the Alven velocity is given by 
VA = J B/(4vrp) 1 / 2 = 4.2 kms^E/lO mG)(n H2 /10 7 cm" 3 )- 1 / 2 = 5.9 kms" 1 . The ratio of the 
turbulent to magnetic energy can be evaluated with /3 tur b = 3(<5t>i os /t> a) 2 = 0.31, implying 
that the magnetic energy dominates over turbulence in this hot core (Girart et al. 2009). 

Unlike the hot core, there are very limited observations addressing the extended struc- 
ture W3(H 2 0)-NW. We are not able to estimate the field strength for B± at the moment. The 
mean polarization angle of (^obs) = 54° ±2° in this extended structure appears quite different 
from the general north-south polarization orientation and implies B± to be at P. A. = —36°. 
Based on its morphology and relative position with W3(OH) and the northern Hn region 
(Figure |2]), we speculate that W3(H 2 0)-NW could be residual matter being pushed by the 
expansion of the northern Hll region and the champagne flow from W3(OH). 

3.3. Comparison with Polarization Measurements on Large Scales 

Our polarization measurements show a decrease of polarization percentage, p, toward 
regions of higher intensity, /, (Figure black dots). To avoid complication caused by the 
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UC Hn region, we exclude polarization measurements in W3(0H). Assuming a power-law 
dependence, we obtain a relationship of logp = (—0.73 ± 0.05) log / + (0.30 ± 0.02) with 
a reduced chi-square of x 2 — 2.2. This "depolarization" phenomenon has been reported 
in many regions, e.g. L1755 (Lazarian et al. 1997), Orion-KL (Rao et al. 1998), OMC-3 
(Matthews, Wilson, & Fiege 2001), W51 (Lai, Girart, & Crutcher 2003; Tang et al. 2009b), 
G5. 89— 0.39 (Tang et al. 2009a). The reduction of polarization may be due to complexity 
of magnetic field geometry (Matthews et al. 2001), decreasing dust alignment efficiency 
(Lazarian et al. 1997), beam smearing over small-scale structures (Rao et al. 1998), and 
integration of varying magnetic domains along the line of sight. 

The SMA polarization measurements nicely fill in the central polarization "hole" in the 
larger scale 850 /im polarization map from the James Clerk Maxwell Telescope (JCMT) of 
Matthews et al. (2009). The three innermost JCMT measurements sit around the edge of 
the 37" SMA field of view. In addition, the JCMT measurements with an angular resolu- 
tion of 20" also show the depolarization effect (Figure [3^; gray triangles). Despite the great 
difference in angular resolutions and map extents between the SMA and JCMT observa- 
tions, the JCMT measurements yield a similar relationship of logp = (—0.72 ± 0.05) log / + 
(1.05 ± 0.03) with x 2 — 1-1. To understand this similarity, we convert the dust intensity 
to column density with assumptions of a gas-to-dust mass ratio of 100, a dust opacity of 
k u = 10 (A/250 yum) _/3 cm 2 g _1 (Hildebrand 1983), and a dust temperature of 20 and 105 K 
for the JCMT and SMA measurements, respectively (Wilson, Gaume, & Johnston 1993; 
Wink et al. 1994). For the central region, we derived (3 = 0.9 by comparing the flux density 
between 2.8 mm (Chen et al. 2006) and 878 /im in maps made with visibilities of projected 
baselines within 8 — 207 kA, the range that the two frequencies have in common. Such small 
is often attributed to grain growth in high-density environment (Miyake & Nakagawa 
1993) and is not valid in the outer region. To convert the JCMT measurements, we adopt 
= 1.6 for active protostellar cores (Rathborne et al. 2010). In Figure |3b, the polarization 
percentage shows a similar dependence with the beam averaged column density between the 
two populations of polarization measurements (Figure [3b). A fit using all the data gives 
logp = (-0.72 ± 0.02) logA^H 2 + (17.7 ± 0.6) with x 1 = 1.4. We suggest that the column 
density may be an important factor when interpreting the depolarization effects caused by 
decreasing alignment efficiency, beam smearing, and integration along line of sight. 

4. Summary 

We present the first interferometric polarization map in the W3(OH) massive star- 
forming region observed with the SMA at 878 /im with an angular resolution of 1"5, corre- 
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sponding to 3 x 10 3 AU. The main findings are summarized as follows: 

1. The magnetic field orientation is well aligned with the synchrotron jet in W3(H20) and 
is close to the plane of sky with B± = 17.0 mG. The total field strength is estimated 
to be 17.1 mG, slightly higher but comparable to the estimate of 10 mG from the 
synchrotron model. The magnetic energy dominates over turbulence in the hot core. 

2. The depolarization effect is discerned in both the SMA and JCMT measurements. We 
find that the polarization percentage has a similar power-law dependence with the 
beam averaged column density for the two populations of polarization measurements. 
The column density may be an important factor to consider when interpreting the 
depolarization effect caused by decreasing alignment efficiency, beam smearing, and 
integration along line of sight. 

3. The extended emission structure in the north-west of W3(H20) has a significant po- 
larization detection that suggests the B± orientation at P. A. = —36°. 

This research is supported by National Science Council of Taiwan through grant NSC 
100-2112-M-007-004-MY2. Vivien Chen would like to thank Dr. B. Matthews for the general 
guideline to convert JCMT SCUPOL flux density. 
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Table 1. Polarization measurements at 878 /xm 
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Table 1 — Continued 



R.A. Offset a 
(arcsec) 


Decl. Offset a 
(arcsec) 


(Jy beam" 1 ) 


P 

(%) 


#obs ± O"6»,obs 

(deg) 


Remarks 


0.78 


-0.59 


0.54 


1.8 ± 0.6 


-70 ± 10 




-0.72 


0.16 


1.81 


0.6 ± 0.2 


-75 ± 9 




0.03 


0.16 


1.98 


0.7 ± 0.2 


-68 ± 7 




0.03 


1.66 


0.38 


2.6 ± 0.9 


1 ± 9 




0.78 


1.66 


0.23 


4.2 ± 1.5 


9 ± 10 





a Position offsets with respect to the map center at (a, 5)(J2000) = 
(2 h 27 m 3.870 s , +61°52'24.60"). 

b The rms of the Stokes I map is 0.027 Jy beam -1 . 



Table 2. Dispersions of Polarization Angles 



Region (0 obB ) 


o0 o bs 0C"6»,obs 


59 


n H2 5v los B ± 


(deg) 


(deg) (deg) 


(deg) 


(10 7 cm- 3 ) (kms" 1 ) (mG) 



W3(H 2 0) 0±2 11.4 6.4 9.4 1.5 1.90 17.0 
W3(H 2 0)-NW 54 ±2 12.0 6.4 10.2 ... 
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Fig. 1. — Polarization map of the W3(OH) region at 878 /im with the total intensity (Stokes 
/) map (contours) overlaid on the linearly polarized intensity, I p , map (colors). Line segments 
give polarization orientation with their lengths proportional to the polarization percentage 
and are plotted where I p > 3uqu and I > 3(7/. The rms value is 07 = 27 mJybeam -1 for 
the Stokes / map and oqu = 3 mJybeam -1 for maps of Stokes Q and U. In the total inten- 
sity map, the western clump is the W3(OH) UC Hn region while the eastern clump is the 
W3(H20) hot core, where triangles mark the two 1.4 mm continuum peaks, A and C, from 
east to west. Two distinct and uniform distributions of polarized emission are observed in 
W3(H 2 0) and W3(H 2 0)-NW. Note that the total intensity in W3(OH) is complicated by the 
free-free continuum of the UC H II region. Additional care is required to derive the polariza- 
tion measurements. Contour levels correspond to (—4, —2, 2, 4, 8, 12, 16, 20, 40, 60, 80) x 07. 
The angular resolution is l'/5. 
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02 h 27 m 04?8 04?7 04?6 04?5 04?4 04?3 04?2 04?1 04?0 

RA (J2000) 

Fig. 2. — The VLA 3.6 cm continuum map (contours; Wilner et al. 1999) overlaid 
on the 878 /im linearly polarized intensity, I p , map (colors). Line segments represent 
the inferred orientation of the projected magnetic fields, B±. Triangles mark the two 
1.4 mm continuum peaks in W3(H 2 0), where the eastern peak, W3(H 2 0)A, is associ- 
ated with the synchrotron jet seen in the VLA 3.6 cm map. Contour levels correspond 
to 30, 60, 120, 180, 240, 300, 360, 500 /iJybeam- 1 . 
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Fig. 3. — (a) Polarization percentage, p, vs. total intensity, I. Black dots are polarization 
measurements with the SMA whereas the gray triangles are measurements with the JCMT. 
(b) Polarization percentage, p, vs. beam-averaged column density, Nn 2 . Despite the great 
difference in angular resolutions and map extents between the SMA and JCMT observations, 
a similar dependence is found with the beam averaged column density. 



